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Abstract 

We address a simple but fundamental issue arising in the study of graphene, as well 
as of other systems that have a crystalline structure with more than one atom per unit 
cell. For these systems, the choice of the tight-binding basis is not unique. For monolayer 
graphene two bases are widely used in the literature. While the expectation values of 
operators describing physical quantities should be independent of basis, the form of the 
operators may depend on the basis, especially in the presence of disorder or of an applied 
magnetic field. Using the inappropriate form of certain operators may lead to erroneous 
physical predictions. We discuss the two bases used to describe monolayer graphene, as 
well as the form of the most commonly used operators in the two bases. We repeat our 
analysis for the case of bilayer graphene. 

1 Introduction 

A peculiar characteristic of graphene is the presence of two atoms per unit cell. The solid- 
state theory for such systems necessitates the introduction of multi-dimensional tight-binding 
bases, whose choice is not unique. The expectation values of physically measurable quantities 
are of course independent of basis; however, in practice this is oftentimes not straightforward 
to see. In particular, if the expectation values of certain operators are to be independent of 
basis, their form must be basis-dependent. 

There appears to exist a rather bit of confusion in the literature about the form of various 
operators in the two tight-binding bases most commonly used to describe graphene. The 
operators that are most commonly misidentified are the fc-space Hamiltonian, the density, 
the density of states, and the single-impurity potential. Some of these operators are used to 
describe the effects of impurity scattering in graphene [1, 2, 3, 4, 5, 6]. Using the correct form 
of these operators is essential in correctly computing the density of states in the presence of 
impurities, which is measured in STM experiments [7, 8, 9]. 

Our purpose is to clarify the subtleties associated with the correct form of these operators. 
We present carefully the two bases, and write down the tight-binding Hamiltonian and its 
low energy expansion in first-quantized language. We also describe the corresponding second- 
quantized formalism, and show that the choice of basis is equivalent to choosing the manner 



of taking the Fourier transform of the second-quantized operators. This allows us to write 
down the form of various operators in the two languages. 

For monolayer graphene, one can choose a basis [10] in which only one point per unit 
cell is used as the origin for the Bloch wave-functions. This basis consists of two p z orbital 
wavefunctions centered on the two carbon atoms of the unit cell; these wavefunctions have the 
same phase factor, determined by the position of the "origin" of the unit cell. Alternatively, 
one can use a second basis, in which the positions of the two atoms in the unit cell are 
used as "centers" for Bloch's theorem; hence the second basis also consists of two p z orbital 
wavefunctions centered at the two carbon atoms, but their phase factors (determined by the 
position of the corresponding atom) are different [11, 12]. 

Bilayer graphene on the other hand has four atoms per unit cell. Consequently, there are 
at least two choices of tight-binding basis. We present the canonical form, which is widely used 
in the literature [13, 14], and in which all four p z orbital wavefunctions have different phases 
(given by the positions of the four atoms in the unit cell). We also discuss an alternative basis, 
in which the four wavefunctions have the same phase factor. 

In section 2 we present the two tight-binding bases and the tight-binding Hamiltonian 
for monolayer graphene and its low energy expansion using a first-quantized formalism and 
Bloch's theorem. In section 3 we present the second-quantized formalism. In section 4 and 
section 5 we present the density operator, and the impurity potential respectively. In section 
6 we discuss the case of bilayer graphene and we conclude in section 7. 

2 Lattice considerations 




a) b) 

Figure 1: Hexagonal honeycomb lattice of graphene (a), and its band structure (b). In b) the equal 
energy contours are drawn, and the Brillouin zone is indicated by dashed lines. The Dirac points K 
and K' are marked by arrows, and the reciprocal lattice vectors a{ 2 are also drawn. 

Given the honeycomb hexagonal lattice of graphene with two atoms per unit cell, one 
can use Bloch's theorem to write down the eigenstates of the lattice Hamiltonian. In the 
tight-binding approximation, one searches for eigenfunctions of the Hamiltonian as linear 
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combinations \l/ fc (r) of atomic wave functions. A common representation of this combination 
is 

tf*(r) = cf(k)Vf k (r) + cf(k)Vf k (r) 

where iV is the number of elementary cells, and the functions 0(r) are the wave-functions of 
the p z orbitals of the carbon atoms. As described below, the coefficients cf^ B are chosen such 
that ^(r*) is an eigenstate of the tight-binding Hamiltonian. The vectors Rj = na\ + ma^ 
with j = (n,m) specify the position of one graphene unit cell, with a\ = av / 3x/2 + 3ay/2, 
and c?2 = — av^3x/2 + 3ay/2, where a is the distance between two nearest neighbors. Also, 
Rj are the positions of the A and B atoms respectively. 

For simplicity we took the positions of the unit cells to be given by the positions of the A 
atoms, 

Rf = R 3 . (2) 
In our choice of the coordinate system, the B atoms are located at Rf = Rj + S 3 , where 

— * — * 

the vector S3 = Sab is one of the three vectors connecting an atom A with its three nearest 
neighbors: Si = a\/3x/2 + ay/2, 82 = —a\^3x/2 + ay/2 and 03 = —ay, as depicted in Fig. 1. 
Note that the choice of the origin, as well as of the axes of the coordinate system is arbitrary, 
but once the choice has been made it has to be used consistently in later analysis. 

In this representation of the tight-binding Hamiltonian eigenstates, one first constructs a 
combination of the atomic wave functions within the unit cell, then attaches a phase factor 
to each cell to construct a Bloch function. This is the "textbook procedure" (see for example 
Ashcroft and Mermin Eq. 10.26 [10]). 

In the second representation one writes the Hamiltonian eigenstates as linear combinations 
of two Bloch functions corresponding respectively to the A and B atoms, but with a different 
phase factor attached to each atom A and B. 



^[e^ A 4(^0(f-4) + e^f 4(%(r - i?f )] (3) 



> 

N 

j 

This second representation is used for example in the paper by Wallace on the band structure 
of graphite [11], and in many recent papers on graphene [12]. 

Note that in each representation we have chosen a tight-binding basis {\l/^ fc (r), ty^ k (r)} 

where v = I /II, and ^f /Bk (r) = j= J2j e^<f>{f - R,f B ), while ^ Bk (r) = j= £\ e^*?* 

—* A I ft 

<p(r — Rj ). We can see that the two bases differ by relative phase factors between their 
components. The eigenstates of the tight-binding Hamiltonian are linear combinations of each 
basis wavefunctions. We will show that, while the coefficients of the linear combinations are 
basis-dependent, the eigenfunctions of the tight-binding Hamiltonian are the same in both 
bases. Also, the expectation value of any physical quantity is independent of the basis chosen. 



3 



2.1 Tight-binding Hamiltonian 



The tight-binding Hamiltonian used to describe graphene allows for hopping between nearest 
neighbors (j,A) and (i,B), such that electrons on an atom of the type A/B can hop on the 
three nearest B/A atoms respectively. Thus we can write 

H = -tY,(\<Pf)(<PF\ + h.c.) , (4) 

(ij) 

where \4>f^ B ) is the standard notation for wavefunctions (4>f^ B \r) = 4>{f — R^ B ). The 
eigenequations for the coefficients c (k) and c B (k) in Eqs.(l,3) are straightforwardly obtained 
from evaluating ((f)^ B \H\^) using Eq.(l), where v = I/II and (^*|r) = \E^(r). Thus we 
obtain 

e(jfc) cj(k) = -t (e- itdl + e~ li ^ + l) cf{k) (5) 
e(jfe) c f(k) = -t (e* Sl + e 1 ^ 2 + l) cf{k) 

in the first basis, or 

e(jfe) c fj(k) = -t (e~ iUl + e~ ith + e~ lU ^ cfj(k) (6) 
e(jfe) c fj(k) = -t (e iUl + e iU * + e iU ^ cfj(k) 

in the second basis. Defining 

f T (k) = -t(e- il31 + e~ ild2 + 1) (7) 

f n (k) = -t(e~ iUl + e~ lth + e- iUs ) , (8) 
the Hamiltonian density is written as (v — I or II) 



/;(*) o 



with the eigenvalues 



e(jfe) = ±|//(jfe)| = ±|/j/(jfe)| = ±ty / 3 + 2cos(v / 3A; ;c a) +4cos(v / 3A; ;c a/2)cos(3A;j / a/2) . (10) 

We should note that the eigenvalues of the Hamiltonian (which give the energy dispersion 
of the two bands of graphene) are the same in both bases, as expected. This is because in the 
two representations, the two functions // and fn differ simply by a phase factor: 

///(*) = fi(k)e- iU - = fi(k)e ik y a , (11) 

— * — * 

where 5ab = 83 = —ay is the vector connecting the A and B atoms in a unit cell. 
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Given that fi(k) = |e(fc)|e l9 ( k \ the Hamiltonian density in the first representation can be 
also rewritten as 

(^, (jE ) o J (12) 

with the phase 9i(k) = — axg[/j(fc)]. 

The k dependence of this phase is shown in Figure (2). One can see clearly the two 
inequivalent Brillouin zone corners K and K* '. Each of the two points is equivalent to all the 
points that can be be obtained by translations with the reciprocal lattice vectors. 
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Figure 2: The k dependence the phase 9j(k) is represented by small segments in the two-dimensional 
k space. One sees clearly the two inequivalent BZ corners K and K f having different topologies, and 
being characterized by opposite Berry phases [15]. 



In the second representation, the Hamiltonian carries an inconvenient phase : 



H n (k) = \e(k)\ 







-i0i i (k) 




(13) 



with 6jj(k) = 9j{k) + k ■ 5ab- The k dependence of the phase 6u(k) is shown in Figure (3). 
We can go back and rewrite the eigenf unctions of the tight-binding Hamiltonian in the two 

bases: 

1 



^2N 



e lk - R > 0(f - Rf) ± e- id 'W(f)(r- Rf) 



(14) 



where the ± signs correspond to the eigenf unctions describing the conduction band, and the 
valence band respectively. In the second representation, 



^ k H {r) = -±= [e iln U(r- R j A ) ± e^^e^f (f>(r - R) 



(15) 



Figure 3: The k dependence of the phase 6u(k) carries an inconvenient addition, so that all the six 
K points of the first Brillouin zone appear different. 



Considering that 9u(k) = 9i(k) + k- 5ab and that Sab = Rf — Rf, one checks easily that the 
two representations lead to the same expression of the eigenfunctions ^ k H {f) = ^f(r). 

2.2 Low energy expansions 

As well known, the energy vanishes at the Dirac points, which are at 



at — aZ 



i?i <- i 2 i -* i ->* 
K mn = £ — 3 — + ma i + na 2 . 

Here £ = ± is the valley index (there are two such points for each elementary cell of the 
reciprocal space). Each point is equivalent to all the points in the reciprocal space that have 
the same £ but different (m, n) and that can be be obtained by translations with the reciprocal 
lattice vectors. The £ = ± pair that is chosen most often contains two corners of the first 

— * — * — * — * 

Brillouin zone, K = K^ and K' = Kq Q , as described in Figs. 1,2,3. Note that 
-> f 27r£ 27r£ 



Thus we can expand the Hamiltonian in the first basis around the Dirac points to find 

fRk) = t 



^VS _ ._ _ N i _ ._ _ . 
f-2~9 ' (°i ~~ °2) ^ 2 g ' ^ + a2j 



where k = K^ n + q and £ = ±1 is the valley index. We see that in this basis the expansion 
does not depend on the choice of (m,n). The six corners of the BZ appear thus equivalent 



to either K or K', and can be recovered by a translation of K and K' by various reciprocal 
lattice vectors. Given the above choice of vectors ai and a 2 , one obtains 

fl(k)=v(Zq x -iq y ) (16) 

where v = 3t/(2a). We can now write the low-energy Hamiltonian density in the 4x4 space 
defined by (KA, KB, K'B, K'A) as: 



H(q) 
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(17) 





( 





e -i6(q) 








\ 






e i9(q) 













v\q\ 













_ e -i9(q) 






V 













/ 



which can be expressed in the compact form: 

Ti(q) = vr z <S> {q x o x + q y (Ty), 

where a and r are the usual Pauli spin matrices. Alternatively we can write the low-energy 
Hamiltonian density in Eq. (17) as: 



Hi(q) 



where 9(q) = arctan(g J/ /g :r ). 

In the second basis, the expression of the Hamiltonian is less convenient because it contains 
the phase factor e~ l ^™ n ^ AB = e % ^ mn '^ a where 

K mn ■ &ab = -2n(m + n)/3 

is independent of the valley index £, but depends on the index (m, n). This makes the six 
corners of the BZ appear inequivalent. Thus in basis //, in the 4x4 space defined by 
(K+ n A,K+ n B,K- n B,K- n A), the Hamiltonian density is 

/ /+ \ 





Hn(q) 



f+* 
J i 



\ 








= v 
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ZmniQx + iQy) 






— * 





(18) 
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where the phase factor z mn = e 2ln ( m + n )/ 3 depends on the choice of the vector K^ n in the 
reciprocal space. In the standard choice for the two valley-points (K = Kq and K' 



oo ) 



we have m 
bases. 



n 



and the low-energy expansion of the Hamiltonian is the same in both 



7 



3 Second quantization 



In the second quantized formalism we can define the operators a], bj that correspond to 

creating electrons on the sublattices A and B, at sites Rf and Rf respectively. From Eq.(l) 
we see that the Fourier transform (FT) of the a, and bj operators should depend on the basis, 
such that 



3 

h{k) = Y,^% (19) 

3 

and 

a II {k) = Y j e rk ^a ] 

3 

M*) = £e*^ (20) 

3 

where the sum is taken over all lattice unit cells. The inverse Fourier transform of these 
operators will be: 



a 3 



[ e-^a^k) = [ e~ iln f a n (k) 

JkeBZ JkeBZ 



bj = / e- lh -^b T (k) = / e~ lk - K i b n (k) (21) 

JkeBZ JkeBZ 

where we define Jg eBZ = J BZ Jr^-, and Sbz = 87r 2 /3y / 3. Given the choice for the origin of the 

unit cell, Rj = Rf, we can see easily that ai{k) = an(k), but bi(k) = e tk ' SAB bn(k). Thus the 
change of basis described in the previous section introduces a different momentum-dependent 
phase factor in the definition of the /c-space Fourier-transformed operators. 

In the second quantized formalism we can write the tight-binding Hamiltonian as: 

H = -tJ2( a h + h-c.) (22) 

<«> 

where t is the nearest neighbor hoping amplitude, and (ij) denotes summing over the nearest 
neighbors. In momentum space the tight-binding Hamiltonian becomes: 

H = [ [aUk)b v (k)f v (k) + h.c] (23) 

JkeBZ 

where v = I /II, and the / functions are defined in Eqs.(7,8) in the previous section. We can 
see that, exactly like in the first-quantized formalism, the form of the Hamiltonian is unique 
in real space, but it depends on the basis in momentum space. 
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4 The density and density of states operators 

It is quite interesting to keep track consistently of the correct form of a few other operators 
in the two bases. We first focus on the local charge-density operator: 

p(r) = YiW- R>h + *<?- Rf) b H ( 24 ) 

3 

In the absence of disorder, the density will be independent of position. However, if impurities 
are present the density will fluctuate, and it is useful to define its Fourier transform: 

p(q) = j d 2 re^p{r) = ^ '"' ^ <>/>., + e * R ' b % ( 25 ) 

J 3 3 

whose expectation value can be related to the results of FTSTS measurements [7, 8, 9]. 

In basis /, the FT of the charge density becomes (see the Appendix for the complete 
derivation): 

p(q) = [ [a\(k) ai (k + g) + e^ B 6}(fc)6j(fc + q)} (26) 

Similarly we can redo the analysis in the basis //: 

P(Q) = i [al^auik + q) + 6} J (fc)6 JJ (fc + q)} (27) 
JkeBZ 

Note that for systems that conserve momentum (translationally invariant), as in the ab- 
sence of disorder and magnetic fields, only the q = term is non-zero. Furthermore, if one 
is interested in average quantities, only the q = term is relevant. Hence, in these cases the 
operators have the same form in the two bases. 

Another operator of interest is the local density of states (LDOS), given by the number of 
electrons of energy a; at a given position. Its integral over uj gives the total density described 
above. The previous formulas can be trivially extended to the local density of states, by taking 
all operators at a specific energy u. 

The expectation values of the density of states operator at various positions on the two 
sublattices and at energy uj are given by: 

(p(Rf,uj)) =11 e-^{a\{k,u3)aj{k + q^)) 

Jq JkeBZ 

= \l e-^tiaUl^anik + lu;)) (28) 

Jq JkeBZ 

and 

(p(Rf,u)) =11 e-^ibKk^Mk + q^)) 

Jq JkeBZ 

= 11 e ^$ 7 (£,u;)M* + ?>)>, (29) 

Jq JkeBZ 
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where fj: eBZ = f BZ jM^- is performed over the first BZ, and /- = / j^z is performed over the 

entire reciprocal space. One can straightforwardly show that if f — Rf only the a) a terms 

contribute, and the tfb terms vanish; conversely, if r — Rf only the tfb terms contribute, and 
the a) a terms vanish. 

Evaluating the density of A and B electrons in the unit cell is also different in the two 
bases. Since in basis / both the A and the B operators are defined at the origin of the unit 
cell, both densities have to be evaluated at this position (which we chose to be the position 
of the A atom, Rf). In the basis //, the density of states is evaluated for each atom at its 
corresponding position (Rf or Rf). 

5 Impurity potential 

We can also write down the form of a delta-function impurity potential. For an impurity 
located on sublattice A, we have 

Va = v A a)a 3 = [ a\(k) ai (k>) 

Jk,k'eBZ 

= f e^'y^\\ l( k)a n (k') , (30) 

Jk,k'eBZ 

while for an impurity on the sublattice B 

Vb = v B b)b 3 = [ e^-^r b \(k) bl (k') 

Jk,k'eBZ 

= I ^'y^b^bu^). (31) 

Jk,k'eBZ 

In the case of a single impurity, it is most convenient to choose the origin of the coordinate 
system such that i?* mp = 0. Thus, in basis / the impurity potential will be independent of 
momentum, regardless of whether the impurity is on the A or on the B site. 

In basis II, Rj imv = 0, and no phase factors will appear when the impurity is on sublattice 
A (at the origin of the coordinate system). However, when the impurity is on sublattice 

~*B ~~ * 

B, R- imp = 5ab, and a momentum- dependent phase factor will appear in the form of the 
impurity potential. We should note that this phase factor comes from choosing the origin of 
the coordinate system on an A atom. Indeed, if one performs a FT of the Friedel oscillations 
generated by an impurity at a B atom while using a coordinate system with the origin at 
a neighboring A atom, one generates a momentum-dependent phase factor in the FT. This 
can be eliminated by changing the origin of the coordinate system from the A atom to the B 
atom, and by carefully tracking the change in the form of the other operators. 

6 Bilayer graphene 

We can generalize the formalism presented in the previous sections to systems with arbitrary 
numbers of electrons per unit cell. Bilayer graphene is made of two coupled graphene mono- 
layers (see Fig. 4), and there are four atoms per unit cell, two for each layer. In real space the 



10 




Figure 4: Bilayer graphene lattice 



tight-binding Hamiltonian can be written as: 




The operators a], b] denote the creation of particles at sites A and B in layer 1, while a], 
b\ denote the creation of particles at sites A and B in layer 2. The sites B in the first layer 
lie on top of the sites A in the second layer, and there is a non-zero t p hopping of electrons 
between them. Also Yl/ij) 12 denotes summing over the nearest neighbors in layers 1 and 2 
respectively; Ylj c denotes summing only over the sites B in the first layer which are on top of 
sites A in the second layer. 

As for monolayer graphene, we can define two types of Fourier transform, consistent with 
two different tight-binding bases. In the first basis one first constructs a combination of 
the atomic wave functions within the unit cell, then attaches a phase factor to each cell 
to construct a Bloch function. The corresponding Fourier transformed operators in second- 
quantized formalism are given by: 



ai(k) 


j 


bj(k) 


j 


a^k) 


j 


h(k) 


= 5>*% 



3 



The vectors Rj = nd\ + ma 2 , with j = (n,m), specify the position of the unit cell of the top 
layer, which we also take to be the origin of the four-atom unit cell of bilayer graphene (see 
Fig. 4). Here we chose the origin of the coordinate system on an A atom in layer 1. 
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In the second basis the positions of the four atoms in the unit cell are used as "centers" 
for Bloch's theorem, and the corresponding Fourier transformed operators are: 



a n (k) 


j 


bn(k) 


3 


a n (k) 


3 


b n (k) 





(34) 

where Rf = Rj and Rf = Rj + 5ab are the positions of the A and B atoms in layer 1, while 

Rf = Rj + Sab and R,f = Rj + 25 ab are the positions of the A and B atoms in layer 2. This 
is the basis that is most often used in the literature to describe bilayer graphene. 
In momentum space the tight-binding Hamiltonian becomes: 

JkeBZ 

+t p a\{k)bi{k) + /i.e.] (35) 

where v = I/II and the fs are the same as the ones defined in Eqs.(7,8) for monolayer 
graphene. 

The density operator is given by 

p{f) = Y^W~ Rf)a]aj + 5(f- Rf)b]bj + 5{r- Rf)a% + 6(r- Rf)b]bj], (36) 

3 

and its Fourier transform is: 

JkeBZ 

where v = I/II, Pi = otj — e^' 5AB , and (3i = e 2t ^' &AB , while in basis II there are no relative 
phase factors, fin = otn = Pn = 1- 

We can also evaluate the density and the density of states at various positions: 

(p(Rf,u;)) =\l e-^(a\(k,u)aj(k + q,u;)) = [ [ e'^t (a^a^k + q,cu)) 



JqJkeBZ JqJkeBZ 

(p(Rf,cv)) =11 e-^ibKk^Mk + q^)) =\l e^f^fi^buik + q^)) 

JqJkeBZ JqJkeBZ 



{p(Rf,u)) = / / e-^icl&ufaik + Zu)) = / / e^t{a\j{k,u)a n {jk + q,u)) 

JqJkeBZ JqJkeBZ 

(38) 



(p(Rf,u)) = // e-^WfauMk + Zu)) = e^(6k(fc,u;)M* + 9>)> 

JqJkeBZ JqJkeBZ 
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where as before Jj: £BZ = f BZ jM^ is performed over the first BZ, while f~ = f is performed 
over the entire reciprocal space. 

Note that (like in the case of monolayer graphene) when working in basis /, all the four 
densities of states are evaluated at the position of the unit cell vector Rj (which we chose to 
be the position of the A atom, Rf). In basis II however, the density of states is evaluated for 
each atom at its corresponding position: R A I B I A I B '. 



7 Conclusions 

We analyzed the two tight-binding bases used to describe monolayer graphene. We showed 
that, while the eigenstates of the tight-binding Hamiltonian, as well as the expectation values 
of physical quantities are independent of basis, the form of certain operators depends on the 
basis. We also showed that the choice of basis is equivalent to the choice of the manner of 
performing the Fourier transforms of second-quantized operators. We wrote down in the two 
languages the Hamiltonian, the density, and the local density of states (LDOS), as well as the 
impurity potential. We also analyzed the case of bilayer graphene and presented two possible 
choices of tight-binding basis, and the form of the aforementioned operators in these bases. 

For the case of the Fourier-transformed density operator, it is important to note that due 
to the (arbitrary) choice of the coordinate system, its expectation value can only be related 
to the FT of the experimental data if this FT is taken using the same coordinate system 
(axes and origin as depicted in Fig.l). If the FT is taken using a different coordinate system, 
a momentum-dependent phase factor is introduced and needs to be accounted for before 
comparing theory and experiment. This inadvertence may lead in some cases to a simple 
rotation of the data, but other more complicated phase factors can also be introduced by a 
mismatch of the origins of the coordinate systems. For the case of a single impurity, it is 
most convenient to use a coordinate system with the origin at the impurity site. However, if 
multiple impurities are present, one needs to keep track of the relative phase factors introduced 
by their spatial distribution. 

We should also comment that our careful tracking of the phase factors generated by the 
change of basis is in general not relevant if one is only interested in uniform properties or 
in spatial averages. However, a careful analysis of the phase factors is crucial if one studies 
systems with disorder, in the presence of an applied magnetic field, and more generally with 
broken translational invariance. 



8 Acknowledgements 

We would like to thank L. Balents, M. Goerbig, F. Piechon, and especially J.-N. Fuchs for 
useful discussions. CB was supported by a Marie Curie Action under the Sixth Framework 
Programme. 



References 

[1] V. V. Cheianov and V. I. Fal'ko, Phys. Rev. Lett. 97, 226801 (2006); V. Cheianov et. al, 
Science 315, 1252 (2007). 



13 



[2] E. Mariani et. al., cond-mat/0702019; C. Bena and S. Kivelson, Phys. Rev. B 72, 125432 

(2005) ; C. Bena, arXiv:0706.4111. 

[3] T. O. Wehling et. al., Phys. Rev. B 75, 125425 (2007); H. P. Dahal, A. V. Balatsky, J.-X. 
Zhu, arXiv:0711.1168. 

[4] N. M. Peres et. al, cond-mat/0705.3040; N. M. Perez, F. Guinea, and A. H. Castro 
Neto, cond-mat/0512091; Phys. Rev. B 73 125411 (2006); M. A. H. Vozmediano, M. P. 
Lopez-Sancho, T. Stauber, F. Guinea, Phys. Rev. B 72, 155121 (2005). 

[5] T. Ando, J. Phys. Soc. Japan 75 074716 (2006); Y. G. Pogorelov, cond-mat/0603327; M. 
I. Katsnelson and A. K. Geim, arXiv:0706.2490. 

[6] Y. V. Skrypnyk and V. Loktev, Phys. Rev. B 73, 241402(R)(2006) and Phys. Rev. B 75, 
245401 (2007). 

[7] P. Mallet et. al, cond-mat/0702406; Phys. Rev. B 76, 041403(R) (2007); G. Rutter et. 
al. Science 317 219 (2007). 

[8] F. Vonau et. al, Phys. Rev. Lett 95, 176803 (2005); F. Vonau et. al. Phys. Rev. B 69, 
081305 (2004); E. Dupont-Ferrier et. al. Europhys. Lett. 72, 430 (2005). 

[9] K. McElroy, et. al, cond-mat/0404005; M. Vershinin, et. al, Science 303, 1995 (2004); 
A. Fang et. al, cond-mat/0404452. 

[10] See for example N. W. Ashcroft, N. D. Mermin, "Solid State Physic" Saunders College 
Publishing (1976). 

[11] P. R. Wallace, Phys. Rev. 71, 622 (1947). 

[12] See for example R. Saito, G. Dresselhaus, M.S. Dresselhaus, Physical properties of carbon 
nanotubes (Imperial College Press, 1999). 

[13] P. Van Mieghem, Rev. Mod. Phys 64,755 (1992); E. McCann, Phys. Rev. B 74, 245426 

(2006) ; E. McCann and V. I. Fal'ko, Phys. Rev. Lett. 96, 086805 (2006). 

[14] J. Nilsson et. al., Phys. Rev. B 73, 214418 (2006); Z.F.Wang et. al. arXivxond- 
mat/06 12483. 

[15] G. P. Mikitik and Yu. V. Sharlai, Phys. Rev. Lett. 82, 2147 (1999). 

9 Appendix 

In basis / the FT of the charge density is written as: 

p(q) = V/ e^' + ^[4(£)a^ (39) 
■ Jk,k'eBZ 

The sum over the lattice unit cells j can be performed to obtain J2rl + Q+Qrl), where 

— * — » 

Qrl is any vector of the reciprocal lattice. However, as the integral over k is constrained 
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to the BZ, not all the terms of the sum contribute to the result, but only those for which 
k! = k + q + Qrl is in the Brillouin zone. For each q and k there is an unique Qrl that 
satisfies this condition. Consequently we have: 

p(q) = f [4(fc)a/(fc + g^+gi ? L) + e^ s fel(fc)6 / (fc + g + Q fiL )] , , (40) 

JkdBZ k+q+Q RL eBZ 

—* —* —* — * — * —* 

However, given the FT definitions in Eq.(19), ai(k + Qrl) = ai(k), and bi(k + Qrl) = bi(k), 
so we have: 

p(q) = [ [4(fc)aj(fc + g) + e^ fl 6j(fc)6j(fc + g)]. (41) 

JkeBZ 

Similarly we can redo the analysis for the second basis: 

=Y,L [e^-^^^^a/z^+e^-^-^ftl/^)^^)] ( 42 ) 
j Jk,k'eBZ 

— * . — * — * — * 

The sum over the sites Rf can be performed to obtain: J2rl^(^ — k' + q + Qrl), where 
Qrl is again any vector of the reciprocal lattice. However, the sum over the sites R? gives 
Y.RL eyi v{-iQRL ■ SAB)S{k -k' + q + Qrl)- Consequently we get: 

P(0) = [ [a\j(k)a H (k + q + Qrl) + e-^^&^M* + q + Q RL )\ , , (43) 

JfceBZ k+q+Q RL eBZ 

—* —* —* — * — * 

From the definitions in Eq.(20) we see that au(k + Qrl) = an(k) and bn(k + Qrl) = 
b n (k)e^ RL ^ AB , and thus 

P(0) = I [a\j(k)a H (k + q)+ + q)} . (44) 

JkeBZ 
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